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The  conditions  under  which  seismic  sensors  offer  opportunities 
for  the  detection  of  acoustic  targets  in  shallow  water  in  the 
1 to  30  Hz  frequency  region  were  examined  using  existing  theo- 
retical and  experimental  knowledge.  Six  aspects  of  the  problem 
were  examined: 

1.  environmental  geometry 

2.  the  physical  environment 

3.  background  noise 

4 . the  geophone  as  a sensor 

5.  seismic  array  design 

6.  measurements  needed  to  improve  knowledge 

Because  of  the  much  higher  attenuation  in  the  bottom  than 
in  the  water,  sound  transmission  in  the  water  is  superior  to 
that  in  the  bottom  when  favorable  sound  velocity  profiles  are 
present  in  the  water.  However,  when  negative  sound  velocity 
gradients  exist  in  the  water,  propagation  through  the  bottom 
offers  the  potential  of  a reliable  propagation  path  that  is 
unaffected  by  the  velocity  gradient  in  the  water.  Thus, 
seismic  sensors  offer  a minimum  range  of  detection.  A 10  to 
20  km  reliable  range  is  estimated  for  seismic  detection  for 
average  bottom  conditions.  Insufficient  experimental  data  exist 
to  confirm  this  estimate;  hence,  several  critical  experiments 
have  been  recommended. 
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Two  issues  arise  in  understanding,  predicting,  and 
sensing  acoustic  signals  propagating  in  the  bottom.  These 
two  are: 

1.  What  is  the  best  method  to  map  the  bottom? 

2.  What  is  the  optimum  way  to  sample  the  propagating 
energy? 

An  analysis  of  the  geophone  as  a reciever  shows  that  its 
inherent  construction  offers  a 4.8  dB  gain  over  an  omnidirec- 
tional hydrophone  in  an  isotropic  noise  field. 

A theoretical  development  was  performed  that  shows  a 
horizontal  geophone  is  superior  to  a vertical  one  for  P-wave 
energy  for  ratios  of  phase  velocities  to  compressional  velo- 
city of  less  than  the  /?. 

Three  short  range  data  sets  were  found  that  show  that  the 
signal  to  noise  ratio  (SNR)  for  horizontal  geophones  is  larger 
than  that  for  hydrophones  located  nearby. 

Quantitative  data  obtained  using  geophones  is  scarce  and 
data  that  compares  the  performance  of  geophones  to  hydrophones 
is  even  more  limited. 

The  attenuation  of  sound  in  the  bottom  will  chiefly 
determine  the  range  to  which  sound  will  propagate  through  the 
bottom  and  yet  have  a useful  amplitude. 
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This  is  the  final  report  of  a study  conducted  in  accord- 
ance with  Task  II  of  ONR  Contract  No.  N00014-74-C-0294 . The  pur- 
pose of  this  study  was  to  identify  conditions  under  which  seismic 
sensors  offer  opportunities  for  detection  of  acoustic  targets 
in  shallow  water,  based  upon  existing  theoretical  and  experi- 
mental knowledge.  Both  signal  transmission  and  noise  were 
considered  in  an  effort  to  determine  whether  sensors  should 
be  embedded  in  the  sea  floor,  be  located  in  the  shallow  water 
environment,  at  some  distance  in  deep  water,  or  on  land. 

In  this  study,  seismic  detection  was  defined  as  the  detec- 
tion of  a signal  in  which  a significant  part  of  the  transmission 
path  is  through  the  ocean  bottom.  Six  aspects  of  seismic 
detection  were  investigated: 

1.  environmental  geometry 

2.  the  shallow  water  environment 

3.  background  noise 

4 . the  geophone  as  a detector 

5.  seismic  array  design 

6.  measurements  needed  to  improve  the  state  of  knowledge. 
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These  are  each  discussed  in  Section  2.  Section  3 discusses 
outstanding  seismic  detection  issues  that  prevent  a definitive 
resolution  of  all  the  problems  investigated  and  outlines  a 
program  for  gathering  needed  information.  Section  4 contains 
a short  summary.  Mathematical  derivations  of  discrete  prop- 
erties of  the  acoustic  sound  field  are  included  in  the  Appendices. 


* 


I 

I 

i; 

i 

r 

i: 

n 


2.  SEISMIC  DETECTION  CONSIDERATIONS 
2.1  Environmental  Geometry 

A generalized  description  of  the  propagation  medium,  or 
the  communication  channel,  together  with  representative  loca- 
tions of  an  acoustic  source  and  receivers  are  illustrated  in 
Figure  1.  The  transmission  medium  and  communications  channel 
analogy  is  intuitively  appealing  since  both  exhibit  the  equiv- 
alent characteristics  of  a non-linear,  memory-less  information 
channel.  The  sources  (targets)  of  interest  are  characterized 
by  (1)  a spectrum  level  of  160  to  180  dB//yPa,  (2)  a 1 to  30 
Hz  frequency  range,  (3)  a source  depth  in  water  of  0 to  400 
meters,  and  (4)  a spectrum  characterized  by  tonals.  The 
typical  locations  for  receivers  that  may  be  used  to  receive 
seismic  energy  are  in  the  bottom,  in  the  water,  or  in  boreholes. 
R-^  and  in  Figure  1 are  located  in  the  water  and  depend  on 
bottom-propagated  energy  to  be  refracted  into  the  water  for 
reception.  R2  and  R^  are  located  at  shallow  depths  in  the 
bottom.  Paths  to  these  receivers  can  be  propagated  either 
through  a bottom  layer  or  from  the  waterborne  energy  near  the 

water-bottom  interface.  Rc  and  R,.  are  located  in  a borehole. 

o 

Rg  is  located  in  the  sediment  and  Rg  in  the  crystalline  base- 
ment. These  two  receivers  are  placed  so  as  to  receive  energy 
propagating  in  the  layer  in  which  they  are  located.  Since  the 
source  (target)  is  always  located  in  the  water  portion  of  the 
communication  channel,  the  total  capability  of  this  time-variable 
heterogeneous  channel  needs  to  be  considered.  The  receivers 
should  be  positioned  to  optimally  sample  the  channel  to  obtain 
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the  best  signal  to  noise  ratio. 

When  favorable  velocity  profiles  exist  in  the  water,  trans- 
mission of  energy  through  the  water  is  superior  to  transmission 
through  the  bottom.  When  the  velocity  profile  in  the  water  is 
unfavorable,  propagation  through  the  bottom  (seismic  detection) 
offers  a potential  of  reliable  propagation  that  is  essentially 
unaffected  by  propagation  conditions  in  the  water.  By  combin- 
ing the  capabilities  of  water  and  bottom  propagation,  the  com- 
munication channel  can  guarantee  a minimum  range  of  detection 
via  predominantly  seismic  paths  and  yet  retain  the  added  capa- 
bility of  the  paths  that  are  predominantly  waterborne. 

The  communication  channel,  as  discussed,  exists  in  both 
shallow  and  deep  water.  However,  it  is  in  shallow  water  that 
it  appears  that  the  largest  gain  in  signal  propagation  can  be 
achieved  by  pursuing  the  total  communication  channel  concept. 

2.2  Shallow  Water  Environment 

Shallow  water  investigations  that  recognized  the  interac- 
tion of  the  water  and  bottom  were  begun  prior  to  World  War  II. 
Probably  the  best  known  participants  in  these  investigations 
were  Ewing,  Pekeris,  Press  and  Worzel  of  Columbia  University. 
Their  work  was  summarized  and  interpreted  after  the  War  in 
Geological  Society  of  America  Memoir  27  (Reference  1) , which 
even  today  is  the  baseline  reference  for  shallow  water  propaga- 
tion. Three  applicable  results  from  the  Memoir  are  the  inter- 
pretation of  sound  transmission;  the  Ewing  effect,  which  is 
the  dispersion  observed  in  the  lower  frequency  region;  and  the 
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concept  of  depth  of  penetration  of  acoustic  energy  into  the 
bottom. 

The  guided  transmission  theory  for  shallow  water  was 
adapted  from  electromagnetic  theory.  Of  particular  importance 
is  the  description  of  the  pressure  and  particle  velocity  as  a 
function  of  depth.  The  vertical  distribution  of  pressure  in 
the  water  has  been  confirmed  by  Ferris  (Reference  2)  and  others. 
Figure  2 is  a typical  amplitude  distribution  versus  depth  for  the 
first  mode  at  a frequency  slightly  above  cut-off,  for  typical 
low  density  and  low  velocity  ratios  for  a liquid  bottom.  Of 
particular  interest  is  the  pressure  amplitude  near  the  bottom. 

For  this  particular  case,  the  receiver  location  for  the  best 
signal  strength  is  near  the  sediment-water  interface. 

The  depth  of  penetration  was  defined  by  Pekeris  (Reference 
1)  as  the  depth  at  which  no  further  information  on  the  struc- 
ture of  the  bottom  at  greater  depths  can  be  obtained  from  dis- 
persion data.  As  a quantitative  measure,  he  used  the  section 
of  the  bottom  that  contained  99%  of  the  total  energy  in  the 
bottom.  Another  way  to  look  at  depth  of  penetration  is  that 
all  useful  guided  energy  will  be  above  that  depth.  Since  we 
are  interested  in  locating  potential  seismic  energy  paths,  the 
depth  of  penetration  tells  us  how  deep  we  have  to  consider  the 
sediments  or  rock.  Attenuation  in  the  sediments  may  also  limit 
the  depths  that  need  to  be  considered.  The  concept  of  penetra- 
tion  depth  was  extended  by  Williams  (Reference  3)  in  his  con- 
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cept  of  hidden  depths. 

The  penetration  depth  for  a given  set  of  environmental 
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parameters  and  a given  frequency  is  an  indication  of  the  amount 
of  energy  available  to  be  propagated  through  the  bottom  to 
large  distances.  One  concludes  from  Figure  2 that  approximately 
one-half  the  propagated  energy  will  propagate  in  the  bottom  for 
this  idealized  situation,  neglecting  attenuation.  In  Appendix  A 
it  is  shown  that  the  penetration  depth  is  inversely  proportional 
to  both  the  differences  between  the  velocity  in  the  bottom  and 
the  phase  velocity,  and  the  thickness  of  the  first  layer  in 
wavelengths.  Thus  one  concludes  that  at  frequencies  appreci- 
ably above  cut-off  where  the  phase  velocity  begins  to  approach 
the  first  layer  velocity,  the  amount  of  energy  coupled  into 
the  bottom  will  become  small  for  the  idealized  situation  of 
two  isovelocity  layers.  For  negative  velocity  gradients  in  the 
water,  there  will  always  be  large  amounts  of  energy  available 
for  propagation  in  the  bottom  and  attenuation  will  be  the 
determining  factor. 

If  another  layer  is  added  to  the  ideal  case  so  that  the 
bottom  is  represented  by  a fluid  layer  and  a substrate,  the 
determination  of  the  depth  of  penetration  appears  to  be  much 
more  elusive.  As  implied  in  Appendix  C,  the  form  of  the  depth 
of  penetration  curve  will  be  the  same  as  for  the  two  layer  case. 
For  the  very  low  frequencies  and  the  high  frequencies  the 
problem  can  be  broken  down  into  appropriate  two  layer  problems 
as  described  in  Reference  4 . 

If  the  bottom  has  elastic  properties,  the  relationship  of 
the  shear  velocity  to  the  liquid  velocity  in  the  layer  above 
determines  the  acoustic  behavior  and  the  above  analyses  are 
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indicative.  When  the  shear  velocity  is  lower  that  the  overlying 
liquid  velocity,  the  situation  appears  to  present  a dichotomy. 

On  the  one  hand,  theory  shows  that  guided  transmission  is  at 
best  minimal,  while  on  the  other  hand  experimental  data  from 
shallow  water  experiments  show  that  guided  transmission  is 
possible  with  varying  amounts  of  attenuation.  This  dichotomy 
may  arise  from  a lack  of  environmental  data. 

The  phenomenon  of  dispersion  in  the  water  wave  was  dis- 
covered by  Ewing  (Reference  1).  Observed  dispersion  is  related 

3 

to  the  depth,  density,  and  velocity  ratios  of  the  transmission 
media,.  Pekeris  makes  the  point  that  if  the  dispersion  curves 
from  different  paths  are  indistinguishable,  so  are  the  parameters 
of  the  media.  From  this  statement  one  would  further  conclude 
that  the  measurement  of  dispersion  is  a pragmatic  way  of  describing 
the  acoustic  behavior  of  the  ocean  bottom. 

In  1955  Kornhauser  and  Raney  (Reference  5)  added  the  effect 
of  attenuation  in  the  bottom  to  the  wave  description  of  propa- 
gation as  a function  of  frequency  and  mode  number,  but  the 
attenuation  coefficients  were  not  specified. 

The  next  significant  step  in  understanding  shallow  water 
propagation  as  well  as  propagation  thrciugh  the  bottom  was  made 
by  Tolstoy,  Clay  and  Blaik  (References  4 and  6) . In  their  Fire 
Island  experiment  it  was  concluded  that  bottom  parameters  deter- 
mined by  seismic  refraction  were  inferior  to  acoustic-determined 
parameters  (Reference  7) . The  experiment  showed  the  adequacy 
of  the  wave  model  to  predict  short  range  shallow  water  propagation. 
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In  Figure  3 relative  propagation  loss  contours  for  a 10  and 
a 14  Hz  source  are  shown  (after  Reference  6) . The  receiver 
is  a vertical  geophone  located  in  a borehole  on  the  beach  at 
a depth  of  119  meters.  The  indicated  propagation  loss  is  greater 
than  would  occur  in  spherical  spreading  (20  log  R) . This  result 
will  be  discussed  further  in  a following  section. 

The  effect  of  a mixture  or  a loaded  liquid  has  been  dis- 
cussed by  Officer  (Reference  8) . For  the  water-bottom  interface 
Figure  4a  shows  the  resulting  velocity  of  the  upper  region  of 
a structureless  bottom  to  be  less  than  that  of  the  water.  Under 
such  a condition  the  natural  question  to  ask  is,  where  is  the 
water-bottom  interface  for  guided  transmission?  The  above 
results  assume  an  ideal  situation.  If  a water-filled  sediment 
is  considered,  then  the  sediment  will  exhibit  some  structure?  i 
probably  very  little  at  the  water-bottom  interface,  but 
increasing  with  depth.  Also  water  has  viscosity.  With 
these  added  real  life  parameters,  consider  a high  porosity 
sediment.  The  viscosity  of  the  water  will  cause  a sediment 
velocity  that  is  dependent  on  frequency  as  illustrated  in  Figure  4b. 
A medium  that  exhibits  such  behavior  will  be  highly  dispersive. 

The  mixture  effects  are  undoubtedly  present  in  the  first  bottom 
layer  and  need  to  be  considered  as  well  as  the  structure  in 
the  sediment. 

In  an  experiment  using  a Texas  Tower  as  a receiving  station, 
Barakos  studied  shallow  waller  propagation  as  a function  of  season 
(Reference  9) . Ranges  in  excess  of  100  nautical  miles  were 
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The  detector  depth  is  389  ft.  The  water  increases  in  depth  perpendicular  to 
the  shore,  and  the  depth  to  basement  increases  in  the  170°  direction.  Attenu 
ation  is  given  for  each  of  the  runs  (represented  by  radial  lines)  in  dB 
relative  to  the  amplitude  at  2 miles.  (After  Reference  6) . 
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achieved  using  bomb  shots  as  sources  with  favorable  sound  velocity 
profiles  in  the  water.  By  using  the  dispersion  concept  of  Ewing 
{Reference  1) , Barakos  acoustically  categorized  a large  area  of 
the  continental  shelf  between  New  Jersey  and  Nantucket.  Bucker 
added  to  the  dispersion  concept  by  advancing  the  idea  of  measuring 
attenuation  by  quantitatively  measuring  the  amplitude  in  small 
frequency  increments  for  each  mode  (Reference  10) . 

To  categorize  the  shelf,  Houtz  and  others  (Reference  11) 
developed  a method  for  measuring  interval  velocities  of  the 
sediments  on  the  shelf  using  sonobuoys.  Interface  waves,  in 
contrast  to  body  waves,  are  a possible  way  to  obtain  acoustic 
information  transmission  over  long  distances.  This  possibility 
is  an  extension  of  teleseismic  observations  at  low  frequencies 
through  the  3 to  30  Hz  region  of  interest.  Two  interfaces  are  of 
special  interest,  the  water-bottom,  and  the  sedimentary  rock- 
crystalline  basement  interfaces. 

Theoretically  there  are  no  limitations  in  the  generation 
of  interface  waves  at  a liquid-solid  interface.  McLeroy  (Ref- 
erence 12)  observed  an  interface  wave  during  his  sea  bottom 
elastic  wave  experiments.  Although  the  wave  did  not  seem  to 
fit  exactly  any  of  the  theoretical  surface  waves,  it  was  inter- 
preted as  a Stoneley  wave.  The  wave  consisted  mainly  of  energy 
in  the  frequency  range  of  6 to  7 Hz,  a frequency  well  below  cutoff 
for  guided  propagation.  Propagation  loss  as  a function  of  range 
indicates  that  the  Stoneley  wave  has  3 dB  less  loss  between  the 
ranges  of  100  and  2000  ft  than  does  the  waterborne  energy. 
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Hence,  the  observed  interface  wave  may  be  a viable  path  for 
distinct  frequencies  (in  this  case  6-7  Hz) . 

In  another  experiment  Bucker,  et  al.  (Reference  13)  used 
Stonely  waves  to  measure  the  shear  velocity  in  the  bottom. 

The  propagating  signal  was  essentially  sinusoidal  at  20  and 
25  Hz  for  each  of  two  cases.  The  extent  of  the  range  measurement 
was  very  limited,  only  47  ft.  The  propagation  loss  was  consid- 
ered to  consist  of  a cylindrical  spreading  loss  and  an  attenuation 
factor  of  0.06  and  0.1  dB/ft  for  each  of  two  cases.  If  the 
attenuation  factor  is  extended  linearly  with  range  to  10  km 
the  loss  would  be  over  3000  dB.  Either  the  propagation  is  very 
poor  or  short  range  measurements  cannot  be  extended  to  longer 
ranges.  In  any  case  more  experimental  data  is  needed  on  surface 
wave  propagation  at  the  water-bottom  interface.  No  information 
was  found  about  a source  in  shallow  water  generating  interface 
waves  between  the  various  rock  layers. 

In  all  the  shallow  water  propagation  studies  of  the  past, 
one  theme  is  highlighted  - an  adequate  acoustic  description 
of  the  bottom  is  needed.  Tolstoy  (Reference  7) , through  a 
sensitivity  analysis,  showed  the  importance  of  accurately  knowing 
the  acoustic  impedance  and  layer  thickness.  An  adequate  descrip- 
tion of  the  bottom  is  one  of  the  prerequisites  for  being  able 
to  predict  acoustic  propagation  in  the  bottom  with  confidence 
and  accuracy.  Another  bottom  parameter  that  has  a controlling 
effect  on  propagation  is  attenuation. 
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Attenuation  in  the  sediments  has  been  of  concern  since 
the  early  1950s.  Tolstoy  developed  empirical  attenuation  coef- 
ficients based  on  his  Fire  Island  data  following  the  method 

of  Kornhauser  and  Raney  (Reference  5) . Nominal  values  for  the 

-4 

exponential  decay  were  5.0x10  /m  at  10  and  14  Hz  (Reference  4). 
Hamilton  and  others  began  measuring  attenuation  and  relating  it  to 
the  porosity  of  sediments.  In  a recent  study  (Reference  14)  he 
summarized  the  attenuation  of  both  compressional  and  shear  waves; 
the  results  are  shown  in  Table  1.  Attenuation  is  in  terms  of 
a logarithmic  decrement  per  wavelength  (A) , which  is  the  natural 
log  of  the  ratio  of  the  amplitude  of  a wave  at  successive  maxima. 
Therefore  attenuation  (a)  as  a function  of  range  can  be  described 
as  a(R)  = e where  R is  range  and  A is  wavelength.  For 

compressional  waves,  sand  and  fine  sand  appear  to  compose  one 
attenuation  group  and  the  silt  clays  and  clayey  silts  another. 

The  difference  in  attenuation  between  these  two  groupings  is  a 
factor  of  six.  Figure  5 shows  attenuation  versus  range  for  three 
values  of  logarithmic  decrement,  0.1,  0.017  and  0.005,  for  a 
frequency  of  3 Hz  and  a wavelength  of  600  m.  The  graph  dramat- 
ically shows  the  high  attenuation  expected  for  a logarithmic 
decrement  of  0.1,  which  is  representative  of  sand.  The  curve 
for  a decrement  of  0.005  is  representative  of  crystalline  rocks. 

Since  there  are  no  known  long  range  measurements  of  prop- 
agation through  these  types  of  sediments,  some  conjectures  will 
be  made  as  to  what  one  might  expect  from  a source  in  the  water. 

A simple  propagation  loss  (PL)  equation  will  be  defined: 
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Table  1 


LOGARITHMIC  DECREMENTS  OF 
ACOUSTIC  ENERGY  IN  SEDIMENTS 


Compress ional 

Waves 

Sheer  Waves 

Low 

High 

Average 

Average 

Sand 

0.073 

0.130 

0.100 

0.3 

Fine  Sand 

- 

- 

0.103 

- 

Silt  Clays 

0.007 

0.028 

0.017 

0.1 

Clayey  Silt 

- 

- 

0.016 

- 

; 


017 


Range  (km) 


PL  = SRL  = aR  + P/S 

where 


SLR  = spreading  loss 
aR  = attenuation 

P/S  = the  ratio  of  compressional  to  shear  energy 

For  the  spreading  loss,  we  will  use  10  log  R + 60  where  R is  in 
kilometers  and  the  60  dB  for  the  first  kilometer  assumes  spherical 
spreading.  The  compressional  to  shear  energy  ratio  in  the  bottom 
sediments  is  derived  from  the  amplitude  functions  of  the  normal 
modes  (Reference  15) . 


K (2-c2/B2)  -K  (z-H) [ 
p/S  = -0 e n 


-H)  [A-c2/a.2  - /l-c2/B2] 


2 (/ 


i 2 

l-c  /a 


where  K is  the  wave  number  of  the  n*"*1  mode 
n 

c = phase  velocity 

a and  & are  the  compressional  and  shear  velocities  in 
the  bottom 
z = depth  variable 
H = water  depth 

The  magnitude  of  the  exponential  decreases  with  depth  and 
is  part  of  the  formulation  to  satisfy  the  condition  of  a vanishing 
potential  at  large  depth.  For  summarizing,  the  exponential 
term  can  be  dropped.  The  attenuation  can  be  obtained  from  Figure  5. 
Let  us  consider  three  types  of  propagation:  (1)  through  the  water 
with  a favorable  velocity  gradient,  (2)  through  a silt  clay  bottom 
and  (3)  through  a crystalline  bottom. 
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For  frequencies  well  above  cut-off  we  would  expect  the 


following  at  a 10  km  range. 


Water  Path  P 
Silt  Clay  Path  P 
Crystalline  Path  S 
♦Unknown 


Type 


P 

P 


SLR 

aR 

P/S 

Total 

PL 

Signal 

Level 

80 

= 6 

0 

86 

160 

80 

23 

* 

103 

160 

80 

22 

1 

103 

160 

Received 

Signal 

Level 


The  above  tabulation  is  probably  optimistic  but  represents  the 
types  of  differences  for  different  paths  through  the  communication 
channel.  The  water  path  is  primarily  controlled  by  the  water 
velocity  gradient,  and  the  boundaries.  When  the  velocity  gradient 
becomes  negative  the  water  path  is  strongly  attenuated  at  large 
ranges.  The  interesting  part  about  the  two  bottom  paths  is 
that  their  propagation  loss  is  about  equal  but  for  different 
reasons.  The  silty  clay  bottom  is  assumed  to  be  poorly  structured 
so  that  it  primarily  supports  P waves.  This  is  not  an  unreasonable 
assumption;  shallow  water  data  interpretations  have  used  this 
assumption  with  reasonable  success  in  the  past.  The  crystalline 
bottom  supports  both  P and  S waves.  However,  the  P/S  ratio 
of  conversion  is  approximately  0.002  for  c equal  to  6 at  3 Hz. 
Therefore  nearly  all  of  the  energy  in  the  bottom  is  in  the  form 
of  shear  waves.  The  attenuation  of  shear  waves  is  approximately 
3 to  5 times  that  of  compressional  waves  resulting  in  a decrement 
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of  approximately  0.025.  Thus  we  have  the  interesting  situation  where 
shear  waves  can  provide  an  effective  means  of  signal  propagation. 
These  numbers  can  be  considered  optimistic,  and  the  controlling 
factor  is  attenuation,  which  increases  linearly  with  range. 

Therefore,  one  can  make  the  following  observations: 

1.  Under  favorable  conditions  the  water  path  will  yield 
useful  signal  levels  at  ranges  of  100  km. 

2.  Under  unfavorable  water  conditions  the  bottom 

path  will  yield  a reliable  range  estimated  to  be  10 
to  20  km . 

2 . 3 Background  Noise 

There  are  few  measurements  of  seismic  background  noise 
in  shallow  water.  Figure  6 shows  typical  observed  noise  levels 
from  the  deep  basins  in  the  western  Pacific  measured  by  ocean 
bottom  seismometers  (OBS)  (Reference  16) . The  particle  velocity 
levels  are  in  general  agreement  with  noise  measurements  using 
hydrophones  in  the  same  general  area  at  similar  depths.  Latham 
(Reference  17)  indicates  that  the  noise  levels  in  shallow  water 
are  similar  to  those  in  Figure  6 except  in  the  frequency  region 
below  15  Hz.  In  shallow  water  the  steep  rise  in  energy  begins 
several  Hz  higher  in  frequency  than  in  deep  water.  Inferences 
from  the  small  amount  of  data  available  can  be  misleading. 

However,  we  can  conclude  that  the  shallow  water  environment 
is  no  noisier  than  the  deep  sea  environment  except  at  the  low 
frequencies.  At  present  we  do  not  know  the  directionality  of 
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Figure  6.  Typical  Noise  Spectra  At  Ocean  Bottom. 
Reference  16) . 
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Data  from  Site  15  (depth,  5600  m)  and  Site  17  (6500  m) . The 
noise  spectra  and  noise  level  seem  to  vary  with  bottom  depth; 
however  when  the  depth  exceeds  5500  m,  the  noise  level  is 
comparable  to  the  quietest  places  on  land,  for  frequencies 
more  than  about  7 Hz. 


the  noise,  whether  it  is  omnidirectional,  primarily  horizontal 
or  vertical.  It  is  not  known  whether  the  directionality  changes 
as  a function  of  frequency.  The  answer  to  these  questions  can 
strongly  affect  the  type  sensor  that  should  be  used. 

Since  geophones  are  particle  velocity  devices,  they  are 
sensitive  to  all  types  of  media  movement.  Any  water  flow  around 
the  sensor  that  causes  the  slightest  movement  can  increase  the 
background  noise  level  significantly.  Wold  and  Purdy  conducted 
an  experiment  in  Buzzards  Bay  in  which  they  had  one  sensor 
mounted  on  a tripod  that  was  resting  on  the  bottom  and  another 
that  was  on  a probe  that  penetrated  the  bottom  (Reference  18) . 
Figure  7 shows  their  results.  The  sensor  on  the  probe  is  much 
quieter  and  shows  no  response  to  the  tidal  currents.  Even  at 
slack  water,  more  than  a 6 dB  reduction  in  background  noise  ap- 
pears to  be  achieved  by  proper  sensor  implantation. 

2 . 4 Geophone  as  a Detector 

A geophone  is  an  electromechanical  sensor  based  on  a non- 
moving reference  point,  and  a moving  earth  that  causes  a voltage 
to  be  generated.  The  mechanical  part  of  the  system  simulates 
a point  fixed  in  space  by  using  a relatively  large  mass.  A 
coil  of  wire  is  wound  on  the  mass.  Permanent  magnets  create 
a magnetic  field  about  the  coil;  the  magnets  are  attached  to 
the  earth.  If  the  earth  moves,  the  magnets  move,  the  coil 
remains  in  place,  and  an  electrical  signal  is  generated  in  the 
coil  proportional  to  the  rate  of  relative  motion;  hence,  the 
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Background  Noise  Measurement  Comparison 
between  Geophones  Mounted  on  a Tripod  and  on  a 
Probe.  (After  Reference  18). 


often  used  term  velocity  geophone.  The  above  ideal  system  is 
difficult  to  achieve.  The  mass  only  approximates  a fixed  point, 
because  it  is  coupled  to  the  earth  by  both  the  suspending  spring 
and  the  damping  dashpot.  The  result  is  a mechanical  resonance 
at  the  low  end  of  the  response  curve,  followed  by  a flat  response. 
The  mechanical  resonance  sets  the  low  frequency  cut-off  for  the 
sensor. 

A geophone  is  a unidirectional  sensor.  An  axial  guide 
restricts  the  motion  of  the  mass;  this  function  is  accomplished 
by  the  spring.  As  a result  of  this  restricted  motion,  only 

I the  component  of  velocity  parallel  to  the  axis  of  the  geophone 

l 

is  sensed.  As  a result,  the  geophone  has  a sensitivity  pattern 
that  is  directly  proportional  to  the  cosine  of  the  angle  measured 
from  the  sensor's  axis.  Figure  8 shows  the  pattern  of  the  geo- 
phone. The  pattern  is  a dipole  with  no  gain  in  any  direction, 
only  attenuation  off-axis.  For  an  omnidirectional  ambient  noise 
field  and  an  on-axis  signal,  the  signal  to  noise  ratio  for  the 
geophone  output  has  been  calculated  to  be  4.8  dB  higher  than 
for  an  omnidirectional  sensor.  The  improvement  in  noise  per- 
formance  should  be  confirmed  and  studied  with  experimental  data. 
Because  motion  restriction  is  accomplished  by  the  mounting  springs, 
some  cross  axis  coupling  can  be  expected. 

Since  a geophone  is  unidirectional,  the  question  arises  as 
to  which  way  it  should  be  oriented  for  optimum  signal  reception. 
Holmer  reanalyzed  the  basic  equations  of  Pekeris  and  arrived  at 


the  conclusion  that  for  a ratio  of  phase  velocity  to  layer  velocity 
of  /T  or  less  the  maximum  particle  velocity  is  in  the  horizontal 
direction  for  compressional  waves  (Reference  Appendix  B) . 

McLeroy  (Reference  12)  has  also  demonstrated  the  same  conclusion 
in  one  of  his  experiments. 

Since  there  appears  to  be  a preferred  direction  for  signal 
arrival  and  the  geophone  possesses  directionality  capabilities 
in  a small  package,  it  is  feasible  to  have  maximum  response 
in  the  signal  direction  and  yet  have  some  noise  reduction. 

To  accomplish  this,  triaxial  or  multiaxial  geophones  can  be  used. 

By  phasing  the  various  axial  components,  they  can  be  electrically 
steered  in  any  desired  direction.  This  capability  can  be  used 
to  increase  the  SNR  at  the  output  of  the  receiver. 

The  coupling  of  a geophone  to  the  medium  is  an  important 
factor  that  often  may  be  overlooked.  The  oil  companies  often 
take  great  pains  to  match  impedances  when  geophones  are  implanted 
in  boreholes.  A dramatic  example  of  the  results  of  different 
degrees  of  coupling  is  shown  in  Figure  9 (after  Reference  18). 

The  two  curves  result  from  the  same  signal  at  the  same  range, 
but  for  different  mounting  and  coupling.  The  upper  curve  is 
the  result  of  a probe  implanted  in  the  bottom  and  the  lower 
is  the  result  of  the  sensor  mounted  on  a tripod  sitting  on  the 
bottom.  The  trace  from  the  implanted  sensor  is  sharp  and 
distinct  in  character,  while  the  other  is  noisy  and  indistinct. 

In  the  light  of  the  directionality  of  a geophone,  let  us 
reexamine  the  Fire  Island  data  shown  in  Figure  3.  Vertical 
geophones  were  clamped  in  the  well  at  a depth  of  119  meters. 
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At  a range  of  2 to  20  miles  one  expects  the  signal  propagation 
vectors  to  be  mainly  horizontal.  Therefore,  if  horizontal 
geophones  had  been  used  it  would  be  expected  that  the  observed 
propagation  loss  would  have  been  somewhat  less,  probably  as  much 
as  10  dB  less. 


2. 5 Seismic  Array  Design 

The  potential  for  deriving  spatial  gain  is  dependent  on 
the  spatial  and  time  coherence  of  a signal.  The  Large  Aperture 
Seismic  Array  has  demonstrated  that  beamforming  is  feasible 
and  results  in  appreciable  gain. 

Considering  the  bottom  as  the  propagation  medium , we  know 
from  facies  changes  associated  with  sedimentation  that  the  prop- 
agation medium  will  be  spatially  dependent.  As  a result  large 
deviations  from  the  plane  wave  assumptions  will  probably  occur, 
and  the  wave  front  curvature  problem  will  need  to  be  effectively 
dealt  with.  In  addition,  the  inhomogeneities  near  the  sensors 
may  cause  serious  deviation,  whereas  the  distant  ones  will 
blend  into  a slowly  varying  simple  medium.  The  inhomogeneity 
effects  near  the  array  can  be  calibrated  in  much  the  same  manner 
that  bottom  arrays  are  presently  calibrated. 

Four  potential  array  designs  will  be  considered:  (1)  the 
distributed  system,  (2)  several  small  aperture  arrays,  (3)  a 
single  large  aperture  array,  and  (4)  borehole  receivers.  Each 
of  these  designs  will  be  considered  in  turn. 
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The  use  of  a distributed  system  is  associated  with  areas 
of  poor  acoustic  propagation.  Although  the  propagation  loss 
for  the  better  propagation  conditions  was  investigated  above, 
there  are  other  water  and  bottom  conditions  in  which  the  prop- 
agation is  much  poorer.  For  these  poorer  conditions,  the  dis- 
tributed system  is  probably  the  ideal  system. 

If  the  distributed  system  is  made  up  of  multiaxial  geophones, 
each  of  these  sensors  can  be  steered  and  an  estimated  spatial 
gain  of  3 dB  may  be  realized.  Thus  each  element  of  the  distributed 
system  is  in  itself  a mini-array. 

Previously  the  reliable  range  for  propagation  through  the 
bottom  was  estimated  at  10  to  20  km  for  silt  clay  type  materials. 
For  sand  bottoms,  the  range  will  be  considerably  shorter.  For  a 
distributed  system,  a sensor  separation  in  the  range  of  10  to 
100  km  can  be  anticipated  depending  on  propagation  conditions. 

Small  aperture  arrays  have  several  outstanding  attributes 
for  use  in  the  shallow  water  environement . Small  aperture  is 
used  here  to  limit  array  size  to  the  range  of  1 to  10  wavelengths. 

With  a small  aperture  the  wave  front  distortion  problem 
should  be  minimized,  thus  normal  beamforming  will  yield  near 
theoretical  gain;  or  if  special  processing  is  required,  it  can 
be  less  complex  than  for  the  large  aperture  array. 

A small  aperture  array  has  greater  spatial  gain  per  installed 
sensor,  since  the  increase  in  array  gain  in  dB  with  increasing 
N is  proportional  to  jjj  . 
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The  use  of  geophones  as  sensors  also  implies  special  pro- 
cedures for  coupling  each  element  to  the  bottom;  therefore, 
the  larger  part  of  array  installation  cost  will  be  the  per- 
element  cost.  The  cost  per  element  will  become  a much  larger 
factor  than  with  hydrophone  arrays. 

Another  aspect  is  that  two  or  more  small  aperture  arrays 
can  be  used  for  multiple  array  processing.  The  type  of  proc- 
essing requires  coherence  between  the  arrays.  The  ambiguity 
plane  processing  gives  high  gain  for  low  resolution  systems. 

Also,  the  multiple  array  processing  yields  excellent  localization 
capability. 

Even  with  high  attenuation  the  use  of  small  aperture  arrays 
in  a multiple  array  processing  scheme  will  have  greater  range 
capability  than  a simple  large  aperture  array,  will  concentrate 
sensors  in  few  locations  as  compared  to  a distributed  system, 
and  will  provide  better  localization  and  classification  infor- 
mation than  the  other  concepts. 

A large  aperture  array  is  useful  for  good  propagation  con- 
ditions and  for  wide  coherent  domain  media.  In  shallow  water, 
propagation  loss  is  in  the  vicinity  of  10  log  R for  optimal 
environmental  conditions,  which  are  estimated  to  occur  less 
than  25%  of  the  time  and  probably  over  less  than  20%  of  the 
eastern  seaboard  continental  shelf.  Under  such  conditions  the 
large  aperture  array  would  make  only  a limited  contribution 
to  surveillance  efforts.  For  example,  a 10  Hz,  50-element 
array  is  approximately  8 km  in  length  for  a compressional 
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velocity  of  1600  m/s.  With  such  a large  aperture,  sophisticated 
calibration  and  beamforming  will  be  required  to  begin  to  realize 
the  potential  spatial  gain.  The  large  aperture  array  does  not 
appear  to  be  a viable  concept  in  shallow  water. 

Locating  geophones  in  boreholes  is  an  interesting  concept 
for  two  reasons:  (1)  the  receiving  location  could  be  near  the 
shore  or  on  shore,  and  (2)  the  propagation  loss  through  the 
bottom  will  be  minimized.  The  minimization  of  the  propagation 
loss  is  based  on  the  higher  attenuation  coefficients  of  the 
unconsolidated  sediments  overlying  the  crystalline  basement. 

The  borehole  concept  would  minimize  the  distance  the  signal 
traverses  the  sediments  to  the  receiver,  which  is  estimated 

to  decrease  propagation  loss  from  several  dB  to  as  much  as  10  > 

to  15  dB.  However,  there  are  difficulties  with  the  concept. 

In  the  Fire  Island  experiment,  the  borehole  was  found  to  resonate, 
and  cause  high  level  artifacts.  It  may  be  possible  to  use 
these  mechanical  resonant  effects  to  advantage.  Little  is  known 
about  the  interaction  of  a signal  with  a borehole  and  its  pos- 
sibilities in  the  3-30  Hz  frequency  region. 

Since  propagation  in  shallow  water  is  controlled  by  the 
size  of  the  attenuation  coefficient,  a distributed  system 
or  small  aperture  array  appears  to  be  the  most  effective  choice. 

^ The  use  of  boreholes  is  a virgin  subject  and  should  be  investigated. 
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2 . 6 Needed  Measurements 


Since  data  that  demonstrates  the  usefulness  of  seismic 
paths  for  acoustic  propagation  at  long  ranges  in  shallow  water 
is  nearly  nonexistent,  all  types  of  measurements  are  sorely 
needed.  However,  it  would  appear  that  concentrating  efforts 
in  the  more  basic  areas  so  that  a firm  foundation  can  be  estab- 
lished upon  which  to  base  surveillance  concepts  may  be  a viable 
basis  for  selecting  needed  measurements.  Therefore,  recommended 
measurements  should  be  limited  to  four  general  areas,  (1)  background 
noise,  (2)  signal  propagation,  (3)  comparative  SNR,  and  (4) 
bottom  acoustic  parameters. 

Background  noise  measurements  should  be  made  at  a number 
of  locations  to  quantify  noise  directionality  and  amplitude  ver- 
sus depth  in  the  communication  link,  i.e.,  from  the  water  surface 
to  some  depth  in  the  bottom.  Comparisons  of  levels  between 
pressure  and  particle  velocity  are  needed. 

Numerous  signal  propagation  measurements  have  been  made 
in  shallow  water,  but  most  of  these  focus  on  the  water  path. 

We  need  some  measurements  that  quantify  propagation  in  the  sedi- 
ments, to  confirm  both  the  wave  theory  and  the  attenuation 
constants,  many  of  which  have  been  determined  at  relatively  high 
frequency. 

Comparative  signal  to  noise  ratios  (SNR)  for  hydrophones  ver- 
sus geophones  need  to  be  determined.  At  short  ranges  there  appears 
to  be  some  benefit  from  using  horizontal  geophones.  Again  these 
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measurements  need  to  be  performed  as  a function  of  depth.  One 
of  the  purposes  should  be  to  verify  the  vertical  velocity  po- 
tential predictions  of  normal  mode  theory. 

The  lack  of  bottom  acoustic  parameter  measurements  has  been 
recognized  as  the  limiting  factor  in  making  propagation  pre- 
dictions in  shallow  water  for  a long  time.  The  situation  has 
not  changed.  However,  the  issue  of  how  to  map  the  bottom  should 
be  resolved  before  proceeding  with  any  extensive  mapping  program. 
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DISCUSSION 


3 . 1 Methodology  Issues 

There  are  many  issues  to  be  resolved  before  seismic  prop- 
agation is  fully  understood  and  its  capability  completely  evaluated. 
However  there  appear  to  be  two  issues  that  are  rather  basic: 

(1)  the  acoustic  description  of  the  bottom  and  (2)  the  sensor 
to  medium  coupling. 

As  previously  stated  the  lack  of  applicable  bottom  informa- 
tion makes  acoustic  signal  predictions  risky  except  in  a quali- 
tative sense.  The  cost  of  mapping  the  bottom  is  expensive. 

At  present,  three  methods  of  determining  the  parameters  are  being 
fostered : 

1.  Seismic  refraction  measurements  of  layer  thickness 
and  velocity  with  density  estimates. 

2.  The  dispersion  method,  which  determines  how  the  bottom 
acts  as  a communication  channel  as  a whole;  only  in- 
homogeneities that  affect  acoustics  are  determined. 

3.  The  inference  method,  which  involves  the  determination 
of  the  effect  of  the  fine  sedimentation  structure 

on  acoustics.  Then  the  sedimentary  sequences  of  the 
bottom  are  mapped  in  detail.  From  these  maps  the 
acoustic  properties  that  one  needs  to  make  predictions 
are  inferred. 

The  first  thing  that  needs  to  be  done  is  to  determine  which 
of  the  above  methods  is  the  most  cost-effective  in  conjunction 
with  acoustic  signal  prediction  accuracy. 


How  to  sample  the  communication  link,  i.e.  the  combined 
water  and  bottom  medium,  is  another  way  of  stating  the  issue 
of  sensor  to  medium  coupling.  Where  should  the  sensors  be  placed? 
Is  the  velocity  potential  function  versus  depth  a good  indication 
of  *~rtical  placement?  The  use  of  hydrophones  versus  geophones 
as  sensors  is  at  present  being  investigated.  Without  under- 
standing the  effects  of  medium  sampling,  the  measured  results 
will  probably  lead  to  contradictory  conclusions.  For  this  issue 
there  are  more  questions  than  answers.  A geophone  versus  hydro- 
phone receiver  comparison  needs  to  be  based  on  physical  reality 
rather  than  on  a set  of  measurements  that  can,  at  best,  be 
related  statistically. 

If  we  can  gain  some  insight  into  these  two  issues,  the  under- 
standing of  seismic  detection  will  take  a significant  step  forward 
and  subsequent  issues  of  importance  may  then  be  examined. 

3 . 2 Critical  Experiments 

The  above  discussion  indicates  qualitatively  that  sound 
propagation  through  the  bottom  may  be  of  importance  to  the  detec- 
tion of  ships.  The  analysis  of  the  distribution  of  the  sound 
potential  as  a function  of  depth  indicates  that  sizable  amounts 
of  energy  can  be  propagated  in  the  bottom.  Short  range  data 
and  theoretical  analysis  indicate  that  the  horizontal  geophone 
may  have  a signal-to-noise  ratio  advantage  over  a pressure  sensor 
at  the  ranges  of  interest.  Since  there  is  no  applicable  data  at 
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these  ranges  of  interest,  three  experiments  are  proposed  to  deter- 
mine the  viability  of  seismic  propagation  in  the  overall  ship  detec- 
tion scheme.  The  guidelines  for  the  experiments  are: 

1.  Both  hydrophone  and  triaxial  geophone  measurements 
should  be  made  for  comparison. 

2.  Although  impulsive  sources  may  be  used  to  define 
propagation  loss  as  a function  of  range,  the  data 
needs  to  be  processed  so  that  it  is  applicable  to  CW 
sources  such  as  ships. 

3.  Basic  environmental  data  needs  to  be  collected  during 
the  experiment. 

4.  All  signal  transmission  paths,  water  and  seismic, 
should  be  measured. 

Experiment  1 - Signal  Propagation  Through  A Silty  Clay  Bottom 

Silty  Clay  is  reported  to  have  the  lowest  attenuation  factor 
among  the  various  sediment  types.  Since  the  eastern  continental 
shelf  is  mainly  composed  of  thick  sections  of  sediments,  this 
experiment  will  investigate  a favorable  acoustic  environment 
for  geophone  sensors.  Signal  propagation  runs  up  and  down  the 
slope  as  well  as  parallel  to  the  slope  should  be  made. 

Experiment  2 - Signal  Propagation  Through  A High  Speed 
Crystalline  Rock  Bottom 

The  purpose  of  this  experiment  is  to  define  propagation  loss 
through  a potentially  low-loss  crystalline  bottom  with  thin 
sediment  cover.  This  experiment  will  evaluate  the  use  of  shear 
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waves  as  an  energy  transmitting  mechanism.  In  general,  the  higher 
the  velocity,  the  lower  the  loss  that  is  expected.  For  crystalline 
rock  the  attenuation  factor  per  cycle  is  very  low,  approximately 
0.05  dB  per  cycle.  This  experiment  will  investigate  a potentially 
favorable  environment  for  long  range  acoustic  propagation  in 
the  bottom  using  geophones  for  sensors. 

These  two  proposed  experiments  will  define  and  quantify 
long  range  signal  propagation  through  the  bottom.  Also,  compara- 
tive signal  to  noise  measurements  should  be  made  between  geophones 
and  hydrophones. 
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4. 


SUMMARY  AND  RECOMMENDATIONS 


From  the  investigation  it  is  concluded  that  there  is 
insufficient  data  available  to  determine  the  impact  that 
seismic  propagation  paths  can  have  on  detection  capability 
in  shallow  water.  The  greatest  contribution  that  a bottom 
path  could  make  would  be  to  furnish  a reliable  range  of 
detection  that  is  independent  of  the  velocity  profile  in  the 
water.  Two  worthwhile  and  critical  experiments  that  would 
better  define  the  potential  role  of  seismic  detection  are: 

1.  Signal  propagation  through  a silty  clay  bottom 

2.  Signal  propagation  through  a high  speed  crystalline 
rock  bottom. 

If  these  experiments  indicate  that  seismic  propagation 
is  a viable  concept,  the  issues  of  bottom  description,  coupling, 
and  hydrophone  versus  geophone  advantages  will  need  to  be 
addressed. 
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APPENDIX  A 

DEPTH  OF  PENETRATION  OF 
ACOUSTIC  ENERGY  INTO  THE  BOTTOM 

By 

Louis  A.  Mole 

The  intent  of  this  Appendix  is  to  provide  a brief  intro- 
ductory explanation  of  the  theoretical  treatment  of  the  depth  to 
which  significant  amounts  of  energy  will  be  present  in  the  bottom. 
The  development  will  follow  the  normal  mode  theory  approach  of 
Pekeris  (Reference  1) . The  principal  observed  features  of  the 
data  records  of  sound  received  from  an  explosion  in  shallow 
water  may  be  explained  by  a model  termed  a two-layer  half  space. 

In  the  model  shown  in  Figure  A-l,  the  upper  layer,  the  water,  is 
a liquid  with  density  and  sound  velocity  c^  whereas  the  lower 
layer,  the  bottom,  is  also  assumed  to  be  a liquid  of  density 
p 2 and  sound  velocity  C2«  The  problem  of  interest  is  to  deter- 
mine the  pressure  field  due  to  an  explosion  in  the  water.  For 
long  ranges,  the  disturbance  can  be  assumed  to  be  produced  by 
a point  source  which,  in  the  absence  of  the  surface  and  the 


Z = 0 


z = d 


point  source 


water  p^, 

z = H ' 

bottom  p2,  c2 

Figure  A-l.  Model  for  a Two  Layered  Half  Space 
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bottom,  would  generate  a spherically  symmetrical  wave  whose 
amplitude  would  decrease  as  the  inverse  power  of  the  range. 
However,  the  initially  spherical  wave  suffers  multiple  reflec- 
tions both  at  the  surface  and  the  bottom,  and  therefore,  at 
long  range,  the  number  of  such  reflections  that  must  be  accounted 
for  is  very  large.  Because  the  requirement  is  to  be  able 
to  compute  the  pressure  variation  at  long  range,  for  a known 
yield  of  explosive  detonated  at  a given  depth  where  the  character- 
istics of  the  bottom  are  known,  this  necessitates  that  an  exact 
solution  of  the  wave  equation  be  attempted. 

Following  Pekeris  (Reference  1) , the  wave  equation  for 
the  sound  potential,  <p,  in  the  water  and  in  the  bottom  is 
written  as: 


0 < z < H (water) 


z > H (bottom) 


(A-l) 


(A-2) 


where  the  subscripts  refer  to  the  water  and  the  bottom,  respec- 
tively. The  acoustic  pressure,  p,  is  defined  as  the  temporal 
derivative  of  the  sound  potential.  The  horizontal  and  vertical 
components  of  velocity,  u and  w,  are  defined  as  the  spatial 
derivative  of  the  sound  potential.  These  equations  are  expressed  as: 


p = p 3<f>  , u = - 3<j>  ,w  = -9<j) 
3t  3r  3z 


(A-3 ) 
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Eqs.  A-l  and  A-2  must  be  solved  subject  to  the  following  con- 
straints: 1)  the  pressure  should  vanish  at  the  surface,  2)  at 
the  source,  the  sound  potential  should  exhibit  a functional 
form  of  f (t-r/c^)/r  where  r is  the  distance  from  the  source 
and  f(t)  is  the  time  variation  of  the  pressure  pulse, 
and  3)  the  vertical  component  of  the  velocity  and  pres- 
sure should  be  continuous  across  the  bottom  interface.  Presented 
formally  these  conditions  are  that: 

<{>^  = 0 at  z = 0 (A-4) 

3c^l  = ^ 2 and  pn<{>-,  = p.,<J>0  at  z = H (A-5) 

3z  9z 


The  solution  for  this  problem  is  obtained  in  two  steps 
by  solving  for  the  case  when  the  point  source  is  periodic  of 
circular  frequency,  w,  and  then  generalizing  the  solution  for 
an  arbitrary  pressure  pulse  f(t)  by  a Fourier  synthesis.  Specif- 
ically, for  an  individual  frequency,  the  sound  potential  is 
written  as: 

(J>  = elu>t  $ (r,z,ui)  (A-6 ) 


and  in  general  it  is  written  as: 

<J>  (r,z,u))  g(u>)  dw 


(r , z , t)  = 1 

2¥ 


t/  eiUt 


(A-7 ) 


where  g(ui)  is  the  spectrum  of  the  pressure  pulse,  defined  as: 
g(w)  = [ e-1(i)t  f(t)  dt  (A-8 ) 
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The  formal  solution  of  the  sound  potential  due  to  a 


periodic  point  source  consists  of  three  equations,  one  for  each 
of  the  regions  in  the  two-layered  half  space.  These  three 
equations  are: 

oo 

sin  3^  cos  B-^H-d)  + ibB2  sin  Bj(H-d) ' fOI 

B^  Bj^  cos  B1H  + ibB2  sin  B-^H  °lzld 


*^=2 


/ 


J 0 (kr ) dk 


*2=2 


J j 0 (kr ) kdk 
0 


sin  B^d 


cos  B1(H-z)  + ib32  sin  ^i<H-zMfor 
cos  B1H  + ibB2  sin " B-^H  J d <_z£H 


and 

Z -iB2(z-H) 

ip  = 2b  J Jn(kr)kdk  Sin  (Bld)  e for  z>H  (A 

0 (B-^  cos  B-^H  + ibB2  sin  B-^H) 

At  each  of  the  boundaries,  Z=d  and  z=H,  the  equations  are  con- 
tinuous. In  addition,  the  following  terms  are  also  defined: 

b = q^/p2  (A-12) 

Bn  = yw2/cn  2 - k2  k < (o)/cn)  (A-13) 


= -i  j/k2  - (w2/cn2)  k > (u)/cn>  (A-14 ) 

Formal  evaluation  of  these  terms  can  only  be  profitably 
made  by  transforming  the  path  of  integration  into  the  complex 
k plane.  At  long  range,  the  only  terms  remaining  from  the  inte- 
gration are  the  residues  from  the  integrands,  which  are  usually 
termed  the  normal  modes.  Therefore,  for  long  range,  the 
asymptotic  form  of  the  solution  may  be  written  as: 
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(A-9 ) 


(A-10) 


11) 
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H / “ V »r  n=l 


i [tot-k  r-  ( it/4 ) ] 

e n F ( x ) 0 (d)  <)>  (z)  (A-15) 

n rn 


where 


) (d)  = sin  (x  d/H) 
n n 


(A-16) 


and 


fl*n>  = 


X 


n 


(A-17) 


2 2 

(x  - sin  x cos  x - b sin  x tan  x ) 
n n n n n 


Eq.  A-15  has  two  forms  depending  on  the  depth  factor  z.  This 
depends  on  the  two  terms: 


b0  = 1 


= b 


0 < z < H 

z > H 


and 


> (z)  = sin  (x  z/H)  0 < z < H 
n n ' 


= sin  (x  ) e 
n 


The  xfi  are  roots  of  the  equations: 


, ~2 — : 2 
x = H.  /a)  - k 


V; 


or : 


tan  x 


n 


-1 


bHV  k 2 


- U) 
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(A-18) 
(A-19 ) 


(A-20) 


kn2  - cj2/c22(z-H)  z > H (A-21) 


(A-22) 


(A-23) 


U 


i (wt-k  r) 

The  factor  e n in  Eq.  A-15  gives  a very  obvious  inter- 

pretation of  k namely  that  k = w/c  , where  c denotes  the 
r n J n n n 

t h 

phase  velocity  of  the  n normal  mode.  The  meaning  of  phase 
velocity  is  that,  in  case  of  an  arbitrary  disturbance,  the 
amplitude  of  the  Fourier  spectrum  of  the  disturbance  at  to  is 
propagated  with  the  speed  of  the  phase  velocity.  The  phase 
velocity  of  the  normal  modes  starts  at  the  cutoff  frequency 
with  the  value  of  the  speed  of  sound  in  the  bottom  and  decreases 
continuously  with  increasing  frequency  toward  the  value  of  sound  \ 

velocity  in  the  water.  For  an  arbitrary  disturbance,  the  low 
frequency  components  in  the  spectrum  of  a pressure  pulse  would  j 

get  ahead  of  the  high  frequency  components,  so  that  at  long 

* 

range,  the  received  pulse  would  appear  in  the  form  of  a train 
of  nearly  sinusoidal  waves  in  which  the  period  decreases  with 
time.  The  factor, 

i [ut  - k r - (tt/4)1 

e n sin(x  z/H)  (A-24) 

n 

in  Eq.  A-15  indicates  that  the  normal  modes  can  be  decomposed  into 
two  plane  waves,  one  traveling  obliquely  upward  and  the  other 
downward. 

The  factor  4>n(z)  in  Eq.  A-15  represents  the  amplitude  of 
th 

the  n normal  mode.  When  c^  > c^  the  roots  of  Eqs.  A-22  and  A-23 

are  real,  whereas  when  c0  < c,  the  roots  are  complex  and  imaginary. 

£ 1 -ik  r 

In  the  latter  case,  the  factor  e n implies  horizontal  attenu- 
ation, whereas  in  the  case  of  a fast  bottom,  there  exist  solutions 
which  suffer  no  horizontal  damping.  For  angles  of  incidence 
greater  than  the  critical  angle,  no  power  is  transmitted  into 
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the  bottom.  For  the  first  mode,  the  amplitude  decreases  expo- 
nentially with  depth  below  the  bottom  and  as  the  frequencies 
become  higher,  on  the  order  of  1 kHz,  very  little  of  the 

I 

energy  in  the  first  mode  is  left  in  the  bottom.  One  consequence 
of  this  weakening  of  penetration  of  the  energy  of  the  first 
mode  into  the  bottom,  with  increasing  frequency,  is  that  ob- 
servations (or  dispersion)  of  high  frequencies  can  yield 
little  information  on  the  nature  of  the  bottom. 

As  a measure  of  the  depth  of  penetration  of  the  first 
mode  into  the  bottom,  a layer  at  the  top  of  the  bottom  is 

defined,  which  contains  99%  of  the  total  energy  in  the  bottom. 

I 

The  significance  of  the  depth  of  pentration  is  that  no  infor- 
mation on  the  structure  of  the  bottom  at  greater  depths  can  ^ 

be  obtained  from  dispersion  data.  The  equation  for  the  pene- 
tration depth  may  be  derived  from  Eq.  A-21  in  the  following 
manner. 

An  expression  will  now  be  developed  to  calculate  pene- 
tration depth  explicitly.  The  depth  at  which  99%  of  the  energy  I 

in  the  bottom  is  above  that  particular  depth  is  equivalent  to 
determining  the  depth  at  which  the  pressure  amplitude  is  0.1 
of  the  value  at  the  bottom  (z=H) . If  the  exponential  portion  of 
Eq . A-21  is  equated  to  0.1,  then  the  logarithm  of  the  result  is: 

-ifi2  (z-H)  = -2.3026  (A-25) 


Using  Eq.  A-14,  the  result  becomes: 


I 

! 


! . If 

I fl 


If  a dimensionless  quantity,  y,  relating  frequency,  depth  and 
velocity,  is  defined  as  follows: 

y = ioH  (A-27 ) 

2ttc^ 

then  Eq.  A-26  can  be  rewritten  as: 

Y (|  - 1)  = (.36647)  (A-28 ) 


Eq.  A-28  can  be  used  to  calculate  the  penetration  depth  for 

specific  discrete  values  of  c where  the  constraint  on  c is 

n n 

that  c^  < cn  < c^.  In  order  to  determine  the  appropriate  value 
of  cn,  the  equation  for  xn  must  be  solved,  i.e.  either  Eq.  A-22 
or  A-23. 

If  a large  depth  of  penetration  is  defined  as  a depth 

greater  than  the  water  depth,  then  for  most  velocity  ratios, 

C1 

— < 0.95,  large  depths  of  penetration  will  be  confined  to 

c2 

areas  with  water  depths  of  less  than  2A.  Hence,  the  lower 
frequency  energy  will  couple  deeper  into  the  bottom. 
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RELATION  OF  PARTICLE  VELOCITY  AND  PRESSURE 
IN  A TWO-LAYER  SHALLOW  WATER  AREA 

By 

Curtis  I.  Holmer 

Geophones  are  inherently  directional  devices  which  respond 
to  the  magnitude  of  vector  particle  velocity  along  the  geophone 
axis  associated  with  acoustic  propagation  in  a medium.  As  a 
result  of  this  directional  response,  a significant  question 
with  regard  to  system  performance  analysis  is  that  of  determining 
an  appropriate  placement  of  the  axis  of  maximum  sensitivity 
for  signal  detection.  The  fact  that  no  universal  answer  to  this 
question  has  appeared  in  the  literature  points  to  the  complexity 
of  the  problem.  In  this  appendix,  we  will  derive  some  general 
answers  to  this  question  for  the  particular  case  of  low-fre- 
quency, shallow-water  propagation  in  a two-layer  half-space, 
through  systematic  analysis  of  an  analytic  model.  While  the 
velocity  is  sensitive  to  the  frequency  range  of  interest,  water 
depth,  bottom  properties  and  source  range,  it  will  be  shown  that 
the  relationship  of  velocity  to  pressure,  as  well  as  the  re- 
lationship between  velocity  components  is  easily  statable. 

The  discussions  to  follow  are  based  on  the  analytic  model  of 
a two-layer  half-space  developed  by  Pekeris  (Reference  1) . 

Liquid  Bottom 

Figure  B-l  shows  the  geometry  of  the  problem  of  interest. 

The  model  assumes  a liquid  (layer  1)  of  depth  H,  density  p.. 
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Figure  B-l.  Assumed  Model  for  a Two-Layered  Liquid 
Half-Space 


constant  sound  speed  overlying  a second  liquid  (layer  2) 
of  density  p2  > p^,  constant  sound  c 2 > c^.  The  acoustic 
source  is  assumed  to  be  in  the  upper  liquid,  and  our  interest 


1 

r 


is  in  the  sound  field  which  occurs  in  either  layer.  This 
model  is  so  highly  idealized  that  dependence  on  it  for  highly 
quantitative  information  is  inappropriate,  but  we  believe 
that  it  is  usefu]  for  answering  qualitative  questions  of 
the  type  of  interest  here. 

The  second  (infinite  depth)  liquid  layer  is  widely  used 
to  represent  the  properties  of  unconsolidated  sediments  and 
is  adequate  for  the  purposes  of  this  discussion. 

Numerous  authors  have  solved  this  particular  problem  in  detail 
so  we  quote  the  appropriate  expressions  from  Pekeris  (Reference  1), 

which  is  widely  available,  rather  than  duplicate  that  analysis.  The 
acoustic  field  may  be  described  by  a velocity  potential  which 
may  be  written  as  the  product  of  the  three  functions  as: 

4>(r,z,t)  = 0 (r)  W(z)  X (wt)  ( B— 1 ) 

from  which  we  can  derive  both  pressure  and  velocity  expressions 
according  to: 

p = p 9_4>_,  u = — 3(|>  , w = — 3 4>  (B-2) 

3t  ~lfr  5z 

where  p is  the  acoustic  pressure,  u and  w are  the  horizontal 
and  vertical  components  of  particle  velocity  respectively.  The 
potential  <J>  for  this  problem  has  the  solution: 

(B-3) 
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where  <f>  ^ and  4>  2 may  be  interpreted  as  a sum  of  cylindrically 
spreading  modes  which  represent  the  energy  which  is  "guided 
in  the  first  layer,  and  have  an  exponentially  decreasing  "pene- 
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tration"  into  the  second  layer.  The  terms  <p  ^ and  <p  2 may 

be  interpreted  as  arriving  from  a spherically  spreading  dipole 

near-field.  Since  we  will  be  interested  primarily  in  longer 

range  propagation,  we  show  here  only  those  long  range  asymptotic 

forms  (valid  for  r > H,  where  the  arrow  indicates  asymptotic 

to  for  large  k r) , i.e.: 
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xn  are  roots  of  the  equations: 
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is  the  horizontal  component  of  mode  wave  number 
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K.  = oj/c.  , K_  = u/c_,  y = / 1 - c.  /c. 


, K2  are  wave  numbers  in  the  infinite  fluids  1 and  2,  not 


mode  wave  numbers. 

The  above  expressions , while  completely  describing  the 
phenomena  associated  with  propagation  in  the  two  layer  liquid 
model,  are  too  complex  for  our  purposes.  We  seek  to  identify 
situations  where  trends  may  be  perceived  more  easily.  First, 
we  observe  that  for  some  large  value  of  k^r,  that: 
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Interpretation 

To  improve  the  visibility  of  the  parametric  dependencies 
described  by  Eqs.  B-7  and  B-8,  we  note  that  the  mode  wave  number 
kn  can  be  viewed  as  a vector  component  of  the  fluid  wave  number 
K^.  Replacing  the  ratio  kn/K^  by  sin  0^  where  0^  = sin-1kn/K^  = 
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sin  '*'(c1/cn).  (The  angle  0^  will  be  interpreted  physically  in  a 
following  paragraph.)  We  observe  then  that  the  expression: 
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can  be  replaced  by  cos  0^.  Eqs.  B-7  and  B-8  are  then  rewritten 
as : 
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’n  z>H  P2C2 
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Manipulation  of  the  terms  in  Eq.  B-lOa  using  the  definitions 
following  Eq.  B-4  will  show  that  this  expression  reduces  to 
Eq.  B-lOb  at z = H,  indicating  that  the  boundary  conditions  of 
continuity  of  vertical  velocity  at  the  boundary  are  preserved. 

The  tangent  term  in  the  denominator  of  Eq.  B-lOa  serves  to 
relate  the  mode  shapes  of  pressure  and  vertical  velocity,  i.e. 
that  minima  in  sound  pressure  produce  maxima  in  velocity,  etc. 

The  expressions  in  Eq S . B-9  and  B-10  suggest  a coherent 
interpretation  of  the  velocity  field  in  the  fluid  as  being  com- 
posed of  that  from  traveling  waves,  with  ray  paths  traveling 
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at  the  angles  + 0 from  the  vertical.  Thus  0 represents  the 
mode  angle.  0 equals  the  critical  angle  for  the  bottom  at  the 
frequency  of  mode  cutoff,  and  increases  smoothly  with  increasing 
frequency  to  90°  from  the  vertical  at  very  high  frequency. 

In  the  water  column,  the  component  of  velocity  in  the  hori- 
zontal direction  is  simply  the  vector  component  in  that  direc- 
tion modulated  by  its  proportionality  to  the  pressure.  The 
vertical  component  is  a similarly  modulated  vector  component 
(but  "out-of-phase",  both  instantaneously  as  indicated  by  the  i 
in  the  ratio  and  in  the  sense  of  having  a maximum  at  pressure 
minima  and  vice  versa) . Turning  to  the  lower  fluid,  we  note 
a similar  vector  component  type  interpretation  of  Eqs.  B-9b  / 

and  B-lOb.  The  magnitude  of  the  horizontal  component  of 
velocity  is  merely  the  magnitude  in  the  upper  fluid  at  the 
interface  scaled  by  the  density  ratio  b = (p^/p2>. 

From  the  above  interpretation,  we  can  see  that  for  mode 
angles  in  the  upper  fluid  which  are  greater  than  45°,  that  the 
velocity  component  in  the  horizontal  direction  in  both  fluids 

will  be  greater  than  those  in  the  vertical.  If  the  critical 

o C2 

angle  for  the  bottom  is  greater  than  45u  (i.e.,  — < 1.41)  then 

1 

the  horizontal  component  is  greater  than  the  vertical  for  all  fre- 
quencies, while  for  large  velocity  contrasts,  we  have  for  low 

c < c2 

frequencies  in  a given  mode  (i.e.,  such  that:  1.41  < — ) 

C1  C1 

that  the  vertical  component  is  larger,  but  a transition  (at 
the  frequency  where  cn  = 1.41  c^)  to  larger  components  in  the 
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horizontal  direction.  In  any  case  the  magnitude  of  both 
components  is  modulated  by  the  magnitude  of  the  pressure  in 
the  mode  at  the  depth  of  interest.  As  a fraction  of  the 
maximum  pressure  in  the  mode,  the  magnitude  of  the  pressure 
at  a constant  depth  in  the  second  fluid  decreases  with 
increasing  frequency  (above  mode  cutoff) . Fig.  B-2  presents 
mode  pressure  distributions  for  the  first  mode  at  several  fre- 
quencies, from  which  the  decrease  in  pressure  for  a constant 
position  in  the  second  fluid  with  increasing  frequency  is 
apparent . 

From  the  above  analysis  we  conclude  the  following: 

• The  components  of  velocity  in  an  underlying  fluid 

in  shallow  water  have  magnitudes  which  are  simply 

related  to  the  pressure  in  the  mode  and  the  mode  angle. 

c2 

• For  small  velocity  contrasts  ( — < 1.41),  or  for  mode 

1 

angles  greater  than  45°  from  vertical  with  any 
velocity  contrast,  the  magnitude  of  horizontal  com- 
ponent of  velocity  is  always  greater  than  the  vertical 
component  at  the  same  depth  in  the  underlying  fluid. 

The  significance  of  this  analysis  for  system  analysis  is 
that  a horizontal  geophone  may  be  used  as  a selective  means 
for  viewing  the  acoustic  energy  in  a mode.  This  may  provide 
operational  advantages  in  terms  of  signal  to  noise  ratio  since 
the  environment  may  operate  selectively  to  reduce  the  magnitude 
of  noises  from  other  sources  relative  to  the  preferred  signal. 
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Figure  B-2.  Vertical  Distribution  of  Pressure  Amplitude  in 
the  Fundamental  Mode  of  the  Free  Wave. 

f = frequency,  in  Hz 
c = velocity  of  propagation  of  free 
wave  in  ft/sec 

. minimum  frequency  = 93.3  Hz 

(After  Pekeris,  Reference  1.) 

(The  example  plotted  here  is  for  p.  = p_  but  similarly  shaped 
distributions  may  be  expected  for  unequal  densities.) 
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Solid  Bottom 

The  influence  of  an  elastic  bottom  on  the  propagation  of 
sound  in  a shallow  water  area  is  significant  when  compared  to 
a liquid  model  bottom  because  of  the  introduction  of  a new 
mode  of  energy  propagation,  the  m = 0 mode,  which  reduces  to 
a Rayleigh  wave  in  the  zero  frequency  limit.  A new  interface 
wave  referred  to  as  a Stoneley  wave,  represents  the  high  fre- 
quency limit  of  this  mode. 

To  adequately  discuss  all  relevant  phenomena,  it  is  desir- 
able to  subdivide  the  discussion  into  two  physically  different 
cases,  depending  on  whether  or  not  a thick  sediment  layer  is 

J 

present.  Taking  the  latter  case  first,  we  have  a shallow  water 
layer  directly  overlying  a basement  rock  structure,  with  little 
or  no  intervening  sediment,  which  we  will  refer  to  as  Case  I. 

Case  II  involves  a shallow  water  layer,  overlying  a shallow  to 
deep  sediment  layer  which  further  overlies  a basement  rock 
structure . 

In  Case  I,  due  to  the  lack  of  a significant  pressure  varia- 
tion component  in  the  sound  speed  profile  in  shallow  water,  the 
sound  speed  variations  through  the  column  are  sufficiently  small 
so  that  the  layer  can  be  treated  as  homogeneous  (as  in  the  liquid) 
bottom  case.  The  following  discussion  is  based  on  the  develop- 
ment of  this  model  by  Ewing,  Jardetsky  and  Press  (Reference  15). 
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The  presence  of  the  solid  bottom  introduces,  through 
coupling  via  Poisson's  ratio  contraction,  two  different  wave 
types  in  the  solid.  These  wave  types  are  described  as  volume 
changing  (i.e.  dilatational , longitudinal  or  P waves)  and 
volume  preserving  (shear  or  S waves) . The  P waves  are  character- 
ized by  alternating  strains  along  the  direction  of  wave  prop- 
agation with  a characteristic  speed  of  propagation  in  the 
solid.  The  analysis  for  the  liquid  bottom  above  is  entirely 
applicable  to  this  P wave  type  for  modes  above  the  first  (m=0) , 
since  the  longitudinal  wave  is  the  same  wave  form  in  a liquid 
or  a solid.  In  addition  to  the  previous  analysis,  the  new 
m=0  mode  type  should  be  added.  This  mode  is  characterized  by 
no  low  frequency  cutoff,  but  a definite  high  frequency  cutoff 
and  transition  to  the  Stoneley  wave  type  at  a frequency  where  the 
ivater  column  is  on  the  order  of  2-5  wavelengths  deep.  Above 
this  frequency  the  mode  propagates  in  the  interface  between  the 
water  and  solid  at  a speed  slightly  less  than  the  speed  of  sound 
in  the  water.  It  is  characterized  by  exponentially-attenuating 
amplitude  away  from  the  interface  in  both  media,  and  implies  a 
mode  angle  of  90°  (from  the  vertical) . Beyond  this  mode,  the 
correspondence  to  the  liquid  bottom  case  is  complete. 

The  S wave  on  the  other  hand  exhibits  material  strains  (and 
hence  response  velocities)  that  are  perpendicular  to  the  direc- 
tion of  wave  motion  and  travel  at  a speed  cg,  which  is  signifi- 
cantly slower  than  c^  for  the  same  material.  For  the  S wave 
type  the  inverse  of  the  analysis  for  the  liquid  bottom  applies 
(i.e.  interchange  cos  0^  and  sin  0^  in  the  expressions  for  u and  w) . 
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The  determination  of  relative  magnitude  of  vertical  and 
horizontal  motion  is  more  complex  because  of  the  need  to  con- 
sider the  relative  strengths  of  P and  S wave  excitation  into  the 
solid  medium,  which  depends  on  the  relative  velocities  of  the 
two  wave  types  in  comparison  with  the  phase  speed  of  a mode  in 
the  water  layer  (i.e.  c /c0  and  c /c  ) . For  Case  I (i.e.  homo- 
geneous  basement  formation  under  a homogeneous  liquid) , and 
typical  basement  properties,  the  speed  of  both  wave  types  is 
much  higher  than  c^  so  that  the  above  ratios  are  small  in  compari- 
son with  unity.  In  this  extreme  both  wave  types  are  excited 
approximated  equally,  with  the  result  that  vertical  and  horizontal 

I 

motions  are  approximately  equal.  The  magnitudes  in  each  of  the 
two  directions  are  equal  with  a 90°  phase  difference,  when  the 
deformation  is  induced  by  the  guided  pressure  *ave  in  the  over- 
lying  fluid.  Both  deformations  in  the  solid  decay  with  an 
exponentially-decreasing  amplitude  away  from  the  fluid-solid 
interface  as  in  the  liquid  bottom  case. 

The  above  simple  picture  does  not  hold  for  solid-fluid 
systems  where  the  fluid  sound  speed  is  comparable  in  magnitude 
with  the  shear  wave  speed  in  the  solid.  The  Case  II  introduced 
earlier  meets  this  latter  situation,  since  the  water  and  sediment 
layers  can  he  modeled  as  an  inhomogeneous  liquid  (i.e.  with  density 
and  sound  speed  varying  through  the  fluid) . The  case  is  of  con- 
siderable academic  interest  and  was  studied  in  detail  by  Tolstoy, 

(Reference  4)  but  we  strongly  suspect  that  it  is  of  little  prac- 
tical interest  for  long  range  detection.  This  assertion  is 
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justified  by  the  considersation  that  close  matching  of  the  wave 
speed  in  the  fluid  with  that  in  the  underlying  medium,  which  is 
necessary  for  good  coupling,  implies  significant  energy  loss 
into  this  medium,  resulting  in  rapid  attenuation  with  range, 
so  that  long  range  detection  will  be  prohibited  by  signal-to- 
noise  ratio  considerations. 

To  summarize  the  qualitative  results  of  this  section,  when 
a fluid  overlies  a homogeneous  elastic  half  space,  and  sound 
speeds  in  the  solid  are  high  with  respect  to  the  fluid,  particle 
velocities  in  the  solid  resulting  from  a signal  propagation  in 
the  fluid  are  equal  in  magnitude  in  both  vertical  and  horizontal 
directions . 

The  above  conclusion  might  then  suggest  that  the  geophone 
orientation  does  not  matter,  but  it  is  clear  that  the  response 
of  a geophone  to  noise  sources  in  a direction  away  from  the  source 
is  not  the  same  for  horizintal  vs  vertical  orientation.  In 
particular,  noise  sources  perpendicular  to  the  geophone  axis  in 
the  horizontal  configuration  would  be  discriminated  against  for 
P wave  direction,  while  this  would  not  be  the  case  for  the  ver- 
tical geophone.  Thus  signal  and  noise  considerations  together 
imply  that  the  horizontal  geophone  is  more  appropriate  for  long 
range  detection,  as  in  the  case  of  the  fluid  bottom. 
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AN  ANALYSIS  OF  SOUND  PROPAGATION 
IN  A THREE-LAYERED  LIQUID  HALF  SPACE 

By 

Louis  A.  Mole 

Although  a large  number  of  authors  treat  the  mathematics  of 
the  three  layered  liquid  half  space,  few  provide  a detailed  and 
explicit  analysis  of  all  aspects  of  theory  of  the  normal  mode 
solution  in  three  layers.  The  intent  of  this  appendix  is  to 
develop  explicitly  the  functional  form  of  the  sound  velocity 
potential  equations  for  each  of  the  three  layers.  The  analysis 
for  the  three  layered  liquid  is  analogous  to  that  of  the  two 
layer  case  as  will  become  apparent  in  the  ensuring  paragraphs. 

The  starting  point  for  the  mathematical  analysis  will  be  the 
definition  of  the  problem  in  the  terms  with  which  it  was 
formulated  by  Pekeris  (Ref.  1).  . 

The  model,  shown  in  Fig.  C-l,  consists  of  three  liquids  which 
are  characterized  by  1)  densities  of  p^,  p2  and  p^,  and  2)  sound 
velocities  of  c^,  c2  and  c3,  respectively.  The  thickness  of  the 
first  layer,  near  the  surface,  is  H and  that  of  the  second 
layer  is  h.  A point  source  of  singular  frequency , w,  is  located 
in  the  first  layer  at  a depth  d below  the  surface.  The  problem 
to  be  analyzed  in  this  appendix  is  the  determination  of  the 
sound  velocity  potential  at  any  point  in  the  half-space.  The 
sound  velocity  potential,  <j> , from  which  the  acoustic  pressure 


and  the  horizontal  and  vertical  components  of  the  velocity  may  be 
determined,  satisfies  the  wave  equation: 
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V2<J>  = 


d2<t> 


(C-l) 


3t 


where  c is  the  sound  velocity  in  one  of  the  three  layers.  In 
order  that  the  solution  be  valid  everywhere  in  the  half  space,  the 
sound  velocity  potential  must  be  zero  at  the  surface  (z  = 0)  and 
asymtotically  approach  zero  at  infinity  (z-*<*>)  . In  addition,  at 
the  boundaries,  z=H  and  z=H+h,  the  vertical  component  of  the 
velocity  and  the  pressure  must  be  continuous.  To  satisfy  these 
conditions,  the  solution  must  be  of  the  form: 
itot 


<J>  = e 


JQ(kr)  F(z)  G (k) 


(C-2 ) 


where  k is  an  arbitrary  parameter  introduced  for  convenience 
for  the  contour  integration  in  the  complex  k plane. 

If  Eq.  C-2  is  substituted  into  Eq.  C-l  the  solution  is  just: 
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The  subscript  n refers  to  the  upper,  middle  or  lower 
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layers,  i.e.,  n = 1,2  or  3.  The  subscript  1 represents  the 
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section  of  the  upper  layer  below  the  source  and  1 represents 
the  section  of  the  upper  layer  above  the  source.  The  boundary 
conditions  for  the  model  require  that: 
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dF. 
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dz 
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with  the  understanding  that  the  solution,  Eq.  C-2,  must  be  inte- 
grated with  respect  to  k from  0 to  00 . 

For  each  of  the  layers,  the  solution  of  Eq.  C-3  requires 

that : 
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B 

sin 

B-^z  + C cos 

61Z 

(C-10) 

2 = 

D 

sin 

&2Z  + E cos 

B2z 

(C-ll) 

3 = 

E 

sin 

-iB3z 

e 

(C-12) 

where  A,  B,  C,  D and  E are  arbitrary  constants  which  depend  on 
the  boundary  conditions  and  strength  of  the  source.  If  Eqs. 

C-6  thru  C-8  are  used  with  Eqs.  C-9  thru  C-12  to  eliminate  these 

I 

constants,  the  appropriate  integral  equations  for  <J>^  , <f>^,  <f>2  and 
can  be  obtained.  However,  for  convenience,  the  following 
definitions  must  be  made: 

C-4 


I 

I 

I 


x = 61H 


(C-13 ) 


b = P±/P2 


(C-14 ) 


g = p2/p3 


s = g33  tan  (32h)  - i82 
g&3  + i62  tan(e2h) 


(C-15) 


(C- 16) 


V = 3^  S cos  x + b32  sin  x 


(C-17 ) 


Then  the  formal  solutions  for  the  sound  velocity  potential 


in  the  various  layers  are: 


<t>3  = 2e 


uo 

f J0 


(kr)  kdk 


sin  3-^z 


(S33  cos  3^  (H-d)  + bB2  sin  3^  (H-d) ] , 0<z<d 


(C-18 ) 


♦ x = 2 eJ 


f J° 


(kr)  kdk 


sin  3-^d 


[S3j  cos  3^  (H-z)  + bB2  sin  3^  (H-z) ] , d<z<H 


. . f sin  B,d 

^2  = 2b  e J J0(kr)  kdk  v~  * 


[S  cos  62  (z-H)  - sin  32(z-H)],  H<z<H  + h 


(C-19 ) 


(C-20) 


oo 

■ . f sin  S,d 

♦ 3 = 2 bg  eia)t  / JQ  (kr)  kdk  


[S  cos  82h  - sin  82h]  e 


-iS3 (z-H-h) 


, z>H  + h 


(C-21) 


I 


______ 


Note  that  Eqs.  C-18  and  C-19  are  identical  except  for  an 
interchange  of  the  variables  z and  d. 


The  solutions  of  the  integrals  in  Eqs.  C-18  thru  C-21 
may  be  obtained  by  evaluation  of  the  appropriate  branch  cuts 
in  the  imaginary  k plane  and  evaluation  of  the  residues. 
However,  in  Ref.  C-l  it  is  demonstrated  that  the  results 
of  the  integrations  in  the  complex  k plane  are  zero  and 
the  normal  mode  solution  is  obtained  from  evaluation 


of  the  residues.  Therefore,  consider  the  equation  for  the 
residues  of  Eq.  C-18,  for  the  upper  most  layer,  which  is  just: 


sin  8-^z 


♦i  " (-2ni)  l (knr)  kn  — ^ 

[S  B,  cos  (H-d)  + be2  sin  (H-d)  ] 

' 9k  ’ 


(C-22 ) 


< 


In  Eq.  C-22,  the  Bessel  function  of  first  order  has  been  replaced 
by  the  Hankel  function  of  the  second  kind,  first  order,  and  the 
appropriate  summation  over  the  variable  kn-  The  derivation  of  the 
expression  for  V represents  the  evaluation  of  the  residues  at 
the  roots  of  the  equation  for  V,  i.e.  V = 0.  The  value  of  the 
derivative  is: 


k [ ] 

&2  cos  eiH 


(C-23) 


where  the  factor  [ ] is: 


[ ] 


- b(B22-81 


2 2 2 
) sin  x cos  x - b&2  x ” hB-^  cos  x * 


2 2 2 2 
(Bj^  + b^B2  tan  x) 


22  2 2 22  2 1 
+ ig  61(B3  -B2  ) cos  x ( B^  + b B2  tan  x) 
* 1 

b3  (b3V  - b2  ) 

And  if  Eq.  C-17  is  solved  for  S when  V is  zero,  the  result 
becomes : 


(C-24 ) 


S = "b62  Sin  61H 

B^  cos  B^H 


(C-25) 


Replacing  the  terms  appropriately  in  Eq.  C-22  yields: 


C-7 


I 


- e^wt (-2ni)  J H<2)(k  rl  kn  sin  61Z 

BlF  "n  I I, 


Q1  B2  COS  B]_H 


(C—  26) 


-b  B2  sin  BjH  ^ cqs  (H_d)  + b B0  sin  6n  (H-d) 

n 77ZTZ  n t7  J-  ^ ■*- 


cos  B-^H 


By  straightforward  manipulation,  the  trigonometric  terms  may  be 

reduced  to  a more  concise  set  of  terms,  such  that: 

' » . iwt  v „ (2)  . k sin  B,z  bB0  sin  B.d 

0 = -2tti  e H (k  r)  _n 1 2 1_  . 

1 n u n Bn  cos  B, H 

1 1 ( 


(C-27) 


Bj^  B2  cos  B^h 


n n 


At  this  point,  it  should  be  noted  that  the  variables  z and  d 
are  symmetric  with  each  other  in  Eq.  C-27,  thereby  obviating 
the  need  for  evaluation  of  Eq.  C-19.  Further  collection  of  terms 


results  in: 


2.i  b e1"1  r H.,z,(k  r)  xn  62  3in  (V/H)  sin  lxn  Z/H>  (C-2S) 

— r; — tO  n 


1 H 


For  the  far  field  approximation,  the  Hankel  function  may 


be  expressed  as: 


ho(2,(V> 


i ( tt/4  - knr) 


Tik  r 
n 


(C-29 ) 


Then  the  formal  solution  for  the  sound  velocity  potential  in 
the  farfield  in  the  upper  layer  is  just: 


* 


>n  = 2 nb  ]~2  ” 1 

"»">/«  ah  FF- 

v n 


i (wt-k  r-ir/4)  „ 2 . , „ 

e n xn  82  sln^-~ 


T I 


d\  lx  z\ 
n_ 1 . ( n ) 

H /Sln\  H /0<z<H 


n 


■ ( C—  3 0 ) 


where  [ ] is  defined  by  Eq.  C-24. 

The  solution  for  the  middle  layer  is  handled  in  a similar 
manner.  The  residues  of  Eq.  C-20  may  be  written  as: 


<j>2  = 2b  e^  < t ( — 2 tt i ) 


l »„,2\ 


k r ) k 
n n 


sin  B , d 


(3Y) 

9k 


l1"*  [S  cos  ^2  (z-H) -sinf^  (z-H)  ] 

(C-31 ) 


n 

As  before,  the  derivative  of  V with  respect  to  k is  evaluated  so 
as  to  yield  the  results  of  Eq.  C-23  and  C-24.  Additionally,  as 
before,  the  roots  of  V are  determined  and  the  expression  for  S, 
Eq.  C-25,  may  be  used.  The  result  for  the  middle  layer  is  then: 


<f>0  = -271  i e1,i)tb  l (k  r)  Kn 

£.  'On  —5 r 7 — • 

n V 'n 


(2) 


k^  sin  B-^d 


(C-32 ) 


-b  32  sin  cos  02 (z-H)  - 8^  sin  B2(z-H)  cos 


8^  cos  B^H 


B1  B2  cos  B^H 

The  trigonometric  terms  may  be  simplified  somewhat  to  yield: 


(“)•“:  "o(2)<v 


sin  B^d 
[ ] ■ ' 


(C-33) 


Xn^2^^2  s*n  xn  cos^2^z_H^  + ^1  cos  xn  s;‘-n  &2^Z-H^ 


C-9 





* 





1 


If  the  farfield  expression  for  the  Hankel  function,  Eq.  C-29 
is  employed,  then  the  formal  solution  for  the  sound  velocity 
potential  in  the  middle  layer  is: 


= [ — 
v2  [ H 


uu 

— I — 

Trr  n=l  /k— 


i(wt-k  r-TT/4)  x 6_  sin  (x  d/H) 
n n z n 


(C-34) 


[b  B0  sin  x cos  B0  (z-H)  + g,  cos  x sin  B0  (z-H)]  H<z<H+h 


The  lower  layer  may  also  be  evaluated  in  the  same  manner  as 
the  previous  two  layers.  The  residues  of  Eq.  C-21  for  the  lower 
layer  may  be  written  as: 


>3  = ~2iri  elwt  bg  £ HQ(2)  (knr)kn 

n 


sin  B^d 


9v  \ 
9k  / n 


(C-35) 


;S  cos  B2h  - sin  B2h)  e 


-iB3 (z-H-h) 


And  again,  if  the  results  for  Eqs . C-23,  C-24  and  C-25  are  employed, 


the  expression  becomes : 

= — 2 7r  i e*,',t  bg  J Hn  ^ (k  r)k  sin  B-,d 
3 ^ L 0 n n 1 

n 

-b  B0  sin  B,  H a . a , 

2 1 cos  B?h  - sin  B?h 

Bx  cos  Bx  H 

_____ 

2 

B3  B2  cos  B-^H 


-iB3 (z-H-h) 


(C-36 ) 


If  the  trigonometric  terms  of  the  expression  are  algebraicly 
manipulated,  the  result  may  be  reduced  to: 


C-10 


<}>2  = -2iTi  e 


iwt 


,-n  sin  ( B-,  d) 

bg  ^ H0  (knr)kn  k [ ] P1P2  e 


-iB^ (z-H-h) 


n n 


(C-37 ) 


[-b  B2  sin  B^H  cos  82^  - 8-^  sin  B2h  cos  B-^H] 


And  then  using  the  farfield  expression  for  the  Hankel  function, 
Eq.  C-29,  the  formal  solution  for  the  sound  velocity  potential  in 
the  third  layer  becomes: 


4>3 


(¥)*/  j 


g x B-, 
^ n 2 


-iB3 (z-H-h) 
e 

(C-38) 


[b  B0  sin  x cos  B^h  + 8,  sin  8->h  cos  x ],  z>H+h 
2 n 2 1 2 n 


The  formal  solutions  for  the  three  layered  liquid  half  space 
are  expressed  in  Eq.  C-30,  C-34  and  C-38  for  the  sound  velocity 
potential  in  the  farfield.  Two  additional  procedures  can  be 
employed  to  validate  these  results.  The  first  is  that  the  value 
of  h can  be  set  to  zero,  thereby  reducing  the  problem  to  a two 
layered  liquid  case.  In  this  instance,  the  expressions  become 
the  same,  and  the  expression  for  the  third  layer  becomes  equal  to 
the  expression  for  the  second  layer  of  the  two  layer  case.  The 
second  check  is  that  the  expressions  for  the  three  layers  may  be 
evaluated  at  each  of  the  two  interface  boundaries,  z=H  and 
z=H+h,  to  test  for  continuity.  When  this  is  done,  the  results 
agree  to  within  the  factor  "b"  of  Eq.  C-14  for  layers  one  and 
two  and  to  within  the  factor  "g"  of  Eq.  C-15  for  layers  two  and 
three,  as  would  be  expected. 
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